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Abstract 

We investigate the dynamics of the sine-Gordon solitons perturbed by spatiotempo- 
ral external forces. We prove the existence of internal (shape) modes of sine-Gordon 
solitons when they are in the presence of inhomogeneous space-dependent external 
forces, provided some conditions (for these forces) hold. Additional periodic time- 
dependent forces can sustain oscillations of the soliton width. We show that, in some 
cases, the internal mode even can become unstable, causing the soliton to decay in 
an antisoliton and two solitons. In general, in the presence of spatiotemporal forces 
the soliton behaves as a deformable (non-rigid) object. A soliton moving in an ar- 
ray of inhomogeneities can also present sustained oscillations of its width. There are 
very important phenomena (like the soliton-antisoliton collisions) where the exis- 
tence of internal modes plays a crucial role. We show that, under some conditions, 
the dynamics of the soliton shape modes can be chaotic. A short report of some of 
our results has been published in [J. A. Gonzalez et al., Phys. Rev. E, 65 (2002) 
065601 (R)]. 
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1 Introduction 



The sine-Gordon solitons are very important in physics. They possess cru- 
cial applications in both particle physics and condensed matter theory. For 
instance, in solid state physics, they describe domain walls in ferromagnets, 
dislocations in crystals, charge density waves, fluxons in long Josephson junc- 
tions and Josephson transmission lines, etc. [1,2,3,4,5,6] 

In general, nonintegrable soliton equations (e.g. the ip A equation and the double 
sine-Gordon) may possess internal degrees of freedom, which are crucial in 
many phenomena [7,8,9,10]. A recent discussion of internal modes of solitary 
waves can be found in Ref. [11]. 

In some cases the internal modes can appear due to the discretization of the 
continuum equations [12,13,14]. 

However, it is well-known that the unperturbed ("pure"), continuum sine- 
Gordon equation does not have internal modes. 

A very remarkable question is the following: can external forces create internal 
modes in the sine-Gordon equation? 

Recently there has been a hot debate in the scientific literature about the 
existence of internal modes of sine-Gordon solitons. 

Some authors [15,16,17,18,19,20,21] have claimed that they have found an 
internal quasimode described as a long-lived oscillation of the width of the 
sine-Gordon soliton. 

On the other hand, a very recent and interesting paper is contradicting all 
these reports [22]. 

By considering the response of the soliton to ac forces and initial distortions, 
Quintero et al [22] show that neither intrinsic internal modes nor "quasimodes" 
exist in contrast to previous reports. We should stress that they use only time- 
dependent perturbations in their work. 

In the present paper we study the sine-Gordon equation perturbed by spa- 
tiotemporal external forces: 

4>tt + 70* - 4>xx + sin0 = F(x,t). (1) 

We will show that with some spatially inhomogeneous forces, the internal 
modes can exist for the sine-Gordon equation. 
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The paper is organized as follows: In Section 2 we investigate the existence 
and stability of the internal modes of sine-Gordon solitons in the presence of 
space-dependent external forces. In Section 3 we present the results of numer- 
ical experiments which confirm the theoretical predictions. In Section 4 we 
study the dynamics of the soliton in the presence of a space-dependent force 
that creates a double-well potential for the soliton. Additionally, the soliton 
is perturbed by a driving spatiotemporal force. Section 5 is dedicated to the 
case where the soliton is moving in an array of inhomogeneities. 



2 Spatially inhomogeneous external forces 



We have shown in previous papers [23,24,25,26,27,28,29] that in equations as 
the following: 

<Ptt + l<Pt - <t>xx + sin = F 1 (x) , (2) 



if the force Fi(x) possesses a zero at x* (Fi(x*) = 0), this can be an equilibrium 
position for the soliton. If there is only one zero, this is a stable equilibrium 
position for the soliton if ( j^^ j > 0. For the antisoliton, it is stable if 

Let us assume that Fi(x) is defined as: 

Fi(x) = 2(B 2 - l)smh(Bx)/ cosh 2 (Bx). (3) 



This is a function with a zero at x* = 0. 

We have chosen this function because it possesses the following properties: (i) 
the exact solution for the soliton resting on the equilibrium position can be 
obtained, and (ii) the stability problem of this soliton can be solved exactly. 
The results obtained with this function can be generalized qualitatively to 
other systems topologically equivalent to this one. Besides, real physical sys- 
tems are related to this example [3]. For instance, in a Josephson junction a 
perturbation that can be described by a function of type F(x) = , where 
R(x) is a bell-shaped function, is an Abrikosov vortex lying in the junction's 
plane perpendicular to its local dimension [30]. 

A similar function can describe a local deformation of a Charge Density Wave 
system [31]. 

Usually, function R(x) is taken as the Dirac's 5-function. However, if we wish 
to model a finite-width inhomogeneity, function (3) is a better choice. 
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The exact stationary solution of Eq. (2), with Fi(x) as defined in (3), is the 
following: 

4>k = 4 arctan [exp (Bx)] . (4) 



The stability analysis, which considers small amplitude oscillations around <pk 
~(j)(k, x) = <f> k (x) + f(x)e xt ] , leads to the eigenvalue problem [23,24,25,26,27,28,29]: 



Lf = Tf, 



(5) 



where L = -d 2 x + fl - 2 cosh^ 2 (Bx)] and T = -A 2 - 7A. 



This problem can be solved exactly [32]. The eigenvalues of the discrete spec- 
trum [23,24,25,26,27,28,29] are given by the formula 



T n = B 2 (A + 2An-n 2 ) - 1, 



(6) 



where A(A + 1) = 2/B 2 . 

The integer part of A (i.e. [A]) yields the number of eigenvalues in the discrete 
spectrum (i.e. n < A), which correspond to the soliton modes: this includes 
the translational mode Tq, and the internal or shape modes T n with n > 
[23,24,25,26,27,28,29]. 

All this theoretical investigation produces the following results (note that pa- 
rameter B will be our control parameter): 

For B 2 > 1, the translational mode is stable and there are no internal modes. 

If I < B 2 < 1, then the translational mode is unstable. However, still there 
are no internal modes. 

When I < B 2 < |, apart of the translational mode, there is one internal 
mode! This internal mode is stable. 

In the case that B 2 < | there can appear many other internal modes! The 
exact number is [A] — 1, where A(A + 1) = 2/B 2 . 

For B 2 < A , 2 . , where A* = 5+ ^ , the first internal mode becomes unsta- 
ble! 
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3 Cauchy problem. Numerical experiments 



What happens when we shift the soliton center- of -mass away from the equilib- 
rium position? 



We have the following initial problem: 
<j){x, 0) = 4 arctan {exp [B (x — x )]} , 
Mx, 0) = 0. 



(7) 
(8) 



In the stable case (B 2 > 1) the center-of-mass of the soliton will make damped 
oscillations (for x ^ 0) around the equilibrium point x* = (see Fig. 1). 




Fig. 1. Soliton center of mass dynamics for the stable case (B 2 > 1). The inset 
shows Fi(x) for B 2 > 1. 



In the case that the translational mode is unstable (| < B 2 < 1), the soliton 
will move away indefinitely from the equilibrium position (see Fig. 2). 



Consider the next initial problem: 

(f)(x, 0) = 4 arctan [exp (Bx)] + C sinh(_Ba;) cosh _A (_Ba;), 

M*, o) = o. 



(9) 
(10) 



In this initial problem the initial soliton is deformed. 
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Fig. 2. Soliton center of mass dynamics for the unstable case {B 2 < 1). The inset 
shows F±(x) for B 2 < 1. 

For | < B 2 < | we will observe oscillations of the soliton width (see Fig. 3). 
This is due to the fact that an internal mode (n — 1) has been excited. 
Eventually, due to unavoidable errors in the initial conditions or to energy 
exchange between the internal mode and the translational mode, the soliton 
will move away from the equilibrium position (remember that the equilibrium 
position is unstable for the soliton center-of-mass) . 



CO 1 




Fig. 3. Soliton's width oscillations when the internal mode can be excited and it is 
stable, \ < B 2 <\. 



It is important to note that the instability of the translational mode does not 
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mean instability of the soliton structure. 

We have to say that the frequency of the oscillations observed in the numerical 
simulations coincides with the one obtained theoretically using equation (6). 
The frequency of the width oscillations can be obtained using the equations 
Ul = y/f\, where r\ = B 2 (3A - 1) - 1. 

The most spectacular phenomenon occurs for B 2 < j-jx~TTj> ^* = 5+ 2^ ■ ^ n 
this case, the first internal mode is unstable. If we study the evolution of the 
soliton from the initial conditions (9)-(10) we will observe the destruction of 
the soliton (see Fig. 4). Two solitons move away (in different directions) to 
"infinity" (or to the boundaries of the system and an antisoliton is formed 
in the place of the original soliton remaining there stabilized. In fact, the 
condition B 2 < 1 implies stability for the center-of-mass of an antisoliton. 




Fig. 4. Soliton's destruction when the internal mode is unstable, B 2 < ^(a.+i) ' 
where A* = S±^Z. 

Note that in this situations, the sine-Gordon solitons do not behave as rigid 
objects, which is what is expected from them in general [33]. 



4 Double-well potential for the soliton 



The initial distortions of the width of the soliton will eventually be damped 
due to dissipation. 

Once the internal modes are possible, as in Eqs. (2)-(3) with | < B 2 < |, we 
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need time-dependent external forces to sustain the oscillations of the soliton 
width. 



On the other hand, if we wish the soliton to remain in some spatially localized 
zone, we need stable equilibrium positions for the center-of-mass of the soliton. 

Let us consider the following spatiotemporal perturbation: 

<t>tt + 70< - <t>xx + sin0 = F 2 (x) + F 3 (x,t), (11) 



where 



Fo(x) 



Fi(x) , if — Xi < x < Xi, 

A n 

- D , if x < -xi, (12) 



cosh [B (x + xi)] 

D A if x>x 

cosh [B (x — xx)] 



andF 3 (x, t) = f cos(ujt)g(x), where g(x) = l/cosh 2 (i?(x-|-a;i))+l/cosh 2 (_E(a;— 
*.))■ 

The space-dependent force F 2 (x) creates a double- well potential for the soliton. 

At the same time, in the interval — X\ < x < x±, we have the same force F^x), 
which was sufficient for the existence of the internal mode. 

Actually, other forces can be used to excite the internal mode. In fact, if we 
have a function F(x) that can mimic approximately the behavior of function 
Fi(x) (specially in the interval —x\ < x < x±, where —x\ and X\ are the 
extrema of function F 1 (x)) when B satisfies the condition | < B 2 < |, then 
this function is good for exciting the internal mode. And note that the behavior 
of function Fi(x) in the interval —X\ < x < X\ is a very common behavior for 
a function in an interval where there is a zero and two extrema. 

If in this double-well potential the middle (unstable) equilibrium position is 
such that the condition for the instability of the first internal mode holds, 
then we can observe again the destruction of the soliton (see Fig. 5). 

As before, an antisoliton is created in the position x = 0. However, in this 
case the "new" two solitons do not move away from the inhomogeneities. 
They remain trapped in the wells. 

The time-dependent force F 3 (x, t) will cause the soliton width to oscillate. The 
center-of-mass of the soliton will also oscillate. 
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Fig. 5. Soliton destruction inside the double-well potential (Eq. (11) with 
B = 0.2, D = 0.0, 7 = 0.1, C = -0.1, f = 0.0, E = 0.2). 

In some cases, the soliton center of mass can oscillate, but these oscillations 
will remain confined inside one of the potential wells (see Fig. 6). On the other 
hand, if we slightly increase the amplitude of the driving force, we can make 
the soliton center of mass to jump between the potential wells created by force 
F 2 (x) (see Fig. 7). 




-20 -10 10 



-10 10 -10 10 
X 



10 10 20 



Fig. 6. Soliton shape dynamics inside the right potential well (Eq. (11) with 

B = 0.5, D = 0.2, x = 0.0, v = 0.0, 7 = 0.1, C = 0.0, / = 0.55, u = 0.55, E = 0.99). 

The snapshots correspond to different time instants. 

Although the soliton is not always in the interval —x\ < x < xi, it will return 
to this interval regularly. While the soliton is in this interval, all the conditions 
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Fig. 7. Soliton shape dynamics when the center of 
mass is jumping between the two wells (Eq. (11) with 
B = 0.5, D = 0.4, x = 0.0, v = 0.0, 7 = 0.1, C = 0.0, f = 0.6, u = 1.0 
E = 0.8). 

hold for the internal mode to be excited. 



A very small (further) increase in the amplitude of the driving force will cause 
extraordinary deformations of the soliton shape (see Fig. 8) . Note that in some 
instances, the soliton can be destroyed, but in this case, this phenomenon is 
not permanent, because the soliton can be later reconstructed. This is a very 
dynamical structure with alternating destructions and reconstructions. The 
center of mass of the structure will be jumping between the potential wells. 

In all these situations the initial soliton is not deformed. All the oscillations 
of the soliton width are produced by the time-dependent external forces. 

We will now change the spatio-temporal force used to drive Eq. (11) with the 
following: 

F 3 (x, t) = f cos(u;t)smh(Bx)g(x) (13) 



where g(x) = l/cosh 2 (_B(x + aq)) + l/cosh 2 (_B(a; — oq)). If / is too large 
(see Fig. 9) the soliton will be immediately destroyed and the "three-particle" 
structure will be a permanent state. There will be oscillations in this structure 
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Fig. 8. Soliton profiles for different time instants governed by Eq. (11) 
(B = 0.5, D = 0.2, 7 = 0.1, f = 0.7, w = 0.55, E = 0.7). Note the large 
deformations in soliton shape. 

but it always will be conformed by an antisoliton and two solitons. 




-40 -20 20 



-20 20 -20 20 
X 



-20 20 40 



Fig. 9. Permanent soliton destruction due to the action of the spa- 
tiotemporal force (13) which is associated to the first shape mode 
(B = 0.5, D = 0.1, x = 0.0, v = 0.0, 7 = 0.1, C = -0.1, f Q = 0.4, u = 1.0). 

When /o is not too large, we can have a situation similar to that of Fig. 8 
(see Fig. 10) where there are large deformations of the soliton profiles with 
sporadic destructions of the soliton and eventual reconstructions. 

The spatiotemporal force F 3 (x, t) in Eq. (11) is chosen such that it acts directly 
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-40 -20 20 -20 20 -20 20 -20 20 40 

X 



Fig. 10. The same as in Fig. 9 but with a smaller amplitude of the driving force 
(B = 0.5, D = 0.1, x = 0.0, v = 0.0, 7 = 0.1, C = -0.1, f = 0.2, uj = 1.0). 

on the internal mode. However, we have corroborated that other spatiotem- 
poral forces also can sustain the soliton width oscillations. This includes the 
"ubiquitous" force F 3 (t) = f cos(ut). 

Here we should say that using the force (13) we can produce harmonic oscil- 
lations of the soliton width (see Fig. 11) even when the center of mass of the 
soliton is oscillating inside a stable potential well. This result was obtained 
using the following equation: 

<t>tt + l<Pt - <t>xx + sin0 = Fi(x) + F 3 (x,t), (14) 

where Fi(x) = _Bsinh( Bx) /cosh 2 (Bx), F 3 (x,t) = focos(ut)smh(Bx)/cosh 2 (Bx). 
Note that in general we will need in all cases spatiotemporal perturbations in 
order to produce oscillations of the shape modes. In particular, in this case 
where the equilibrium position is stable, a spatiotemporal force that acts di- 
rectly on the first internal mode is required. 

Now we will return to the system where the soliton is moving in a double-well 
potential created by an external space-dependent force F 2 (12): 

<ptt + 70< - <\>xx + sin0 = F 2 (x) + F 3 (x, t), (15) 

where F 3 (x,t) = f cos(cut)g(x), where g(x) = l/cosh 2 (_E(x+Xi))+l/cosh 2 (_E(a;— 
xi)). In addition, B = 0.5, D = 0.4, 7 = 0.1, f = 0.6, u = 1, and E = 0.8. 
Here the parameters are chosen in such a way that the soliton center of mass 
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Fig. 11. Time series of the soliton width for a one well potential, see Eq. (14) 
(B = 0.5, a = 0.0, f = 0.1, u = 1.0). 



will have a dynamics similar to a Duffing equation [34]: 

d cc doc o . . 

— + + X = /o cosM). 



(16) 



In this case, the soliton center of mass will perform chaotic oscillations jumping 
between the two wells created by force F 2 (see figures 12, 13). However, the 
most interesting result here (considering that our main subject is the internal 
mode) is the fact that we can observe chaotic oscillations of the soliton width 
(see Fig. 14). 



5 Disordered systems 



The propagation of solitons in disordered media has been studied intensively 
in recent years [29,35]. 



Consider the equation 

<Ptt + l<Pt ~ <Pxx + sin0 = F(x). 



(17) 
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Fig. 12. Phase portrait produced using the dynam- 
ics of the soliton center of mass governed by Eq. (15) 
(B = 0.5, D = 0.4, x = 0.0, v = 0.0, 7 = 0.1, C = 0.0, f = 0.6, u = 1.0, 
£ = 0.8). 
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Fig. 13. Poincare map for the dynamics of Fig. 12 
(5 = 0.5, D = 0.4, x = 0.0, w = 0.0, 7 = 0.1, C = 0.0, f = 0.6, u = 1.0 
£ = 0.8). 



where F(x) is defined in such a way that it possesses many zeroes, maxima 
and minima (see Fig. 15). This system describes an array of inhomogeneities. 
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Fig. 14. Poincare map produced using the dynamics of the soliton width (see Eq. 
(15)) (B = 0.5, D = 0.4, x = 0.0, v = 0.0, 7 = 0.1, C = 0.0, f = 0.6, u = 1.0 
E = 0.8). 




Fig. 15. Array of inhomogeneities (see Eq. (18)). 
For our study, we have denned F(x) in the following way: 

e B(x+x n ) _ e 3B(x+x n ) 



F(x) = £ 4 (l - B 2 ) 



n=—q 



^p.B(x+x n ) _|_ iy 



(18) 



where x n = (n + 2) log yy2 + lj / B (n = —q, —q + l---,q — l,q), and q + 2 
is the number of extrema points of F(x). 
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In our array, there is a "superposition" of the "disorder" with a dc component 
which will cause the soliton to move to the right all the time. 

When the soliton is moving over intervals where < 0, the internal mode 
can be excited. In fact, the points Xi where F(xi) = and dF ^ < 0, are 
"barriers" which the soliton can overcome due to its kinetic energy. These 
"collisions" with the barriers will excite the internal modes if in these intervals 
the function F(x) mimics the behavior of F±(x) when | < B 2 < |. 

In Fig. 16 we can see how the width of the soliton will perform sustained 
oscillations during its motion in a disordered medium. 




(0 



Fig. 16. Fourier transform of the dynamics of the soliton width while it is moving 
in the array of inhomogeneities (see Eqs. (17) and (18)) 



6 Discussion and conclusions 

We have shown that in sine-Gordon equations perturbed by inhomogeneous 
(space-dependent) forces F(x), the solitons can possess internal modes. 

Our results are in agreement with previous works [22,33]. In fact, in Eq. (2), 
with Fi(x) as in (3), if we put B 2 = 1, then there are no external forces, and 
from Eq. (6) we obtain that there are no internal modes in that case neither. 

Moreover, having any external inhomogeneous force does not necessarily mean 
having internal modes. 

For instance, if F(x) has a zero which corresponds to a stable equilibrium 
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position for the soliton, even then the internal modes are impossible. This 
explains why it has been so difficult to find sine-Gordon internal modes. 

For the existence of internal modes in the sine-Gordon solitons we need zeroes 
of the function F(x) that corresponds to unstable positions for the center-of- 
mass of the soliton. 

When the soliton center-of-mass is very close to an unstable equilibrium po- 
sition, there is a pair of forces acting in opposite directions on the "body" 
of the soliton. This pair of forces should be sufficiently large to stretch the 
soliton "body" , such that the soliton internal mode can be excited. 

Function F(x) can possess several zeroes corresponding to unstable and stable 
equilibrium positions. For instance, we have studied a force F(x) that creates 
a double-well potential for the soliton. 

Periodic time-dependent forces (along with the spatially inhomogeneous forces) 
can sustain the oscillations of the soliton width. 

A soliton moving in an array of inhomogeneities can also undergo sustained 
oscillations of its width. All this is possible because the internal mode of the 
soliton can exist when it is moving in media where there are inhomogeneous 
space-dependent forces with unstable equilibrium positions. 

Nonetheless, we have discovered another more remarkable phenomena. The 
sine-Gordon internal mode not only can exist for some external forces, but (in 
some situations) it can become unstable. 

If we have an unstable equilibrium position for the soliton center-of-mass and 
the pair of forces acting on the soliton is too large, then the soliton can be 
destroyed. When the soliton is destroyed, it can be transformed into an an- 
tisoliton and two new solitons. The topological charge is conserved. We had 
found this phenomenon before for the 4 -equation [24,25]. However, here for 
the first time we have shown not only that the sine-Gordon soliton internal 
mode can exist, but that it can become unstable and destroy the soliton. This 
is a spectacular manifestation of the fact that the sine-Gordon soliton can 
behave as a deformable (non-rigid) object. 

We should say here that an effective potential for the soliton with stable 
and unstable equilibrium positions can be created also with the following 
perturbation 

<Ptt + 70* - (fixx + sin = H(x) sin 0, (19) 

where the extrema of function H(x) can correspond approximately to the 
equilibrium positions for the soliton (i.e. not the zeroes as in the case of force 
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F(x)). With some appropriate functions H(x), the internal mode of the soliton 
can be excited too. 

In Ref. [3], and the references cited therein, we can find many physical appli- 
cations for such a system. 

It has been recognized that the internal mode is a fundamental concept for 
the nonintegrable soliton models [11]. The internal modes of solitons govern 
the resonant soliton collisions and the soliton- impurity interactions [7,8,9,10]. 

The question about the conditions for creating the soliton's internal mode has 
interested physicists very strongly [11,12,13,14,22]. Several authors [11,12,13,14] 
have shown how the internal mode can appear due to either a small deforma- 
tion of nonlinearity or a weak discreteness. However, this is the first time that 
external perturbations are shown to be able to create an internal mode for the 
sine-Gordon soliton. 

The existence of the internal mode is very relevant to the soliton dynamics 
because it can temporarily store part (or all) of the soliton's kinetic energy, 
which can later be re-stored again in the translational mode. This is important 
in the soliton-antisoliton collisions [7,8,9]. 

The coupling between the translational and the internal modes is the governing 
factor in many contexts, such as the interaction with inhomogeneities, with 
thermal noise, and with external drivings. 

The soliton internal shape modes can contribute to the thermodynamic prop- 
erties of the collective soliton-phonon gas [36,37]. 

Considering all these facts, we believe that the proof of the existence of sine- 
Gordon internal modes in the presence of spatiotemporal perturbations is a 
fundamental result. 
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